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Abstract 

^ \ We perform a calculation of the absolute charged current neutrino-oxygen events 

^— ^ ' rates relevant in the atmospheric neutrino experiments. The inclusive reaction cross- 

pg ■ section is split into exclusive channels, which are classified according to the number of 

CN ! Cerenkov rings they produce. The model includes the effects of residual interaction in a 

RPA scheme with both nucleon-hole and Delta-hole excited states and the effects of [np- 
Q^ ■ nh) excitations {n=2,3). Our result is that although the flavor ratio ^/e remains almost 

unaffected by the nuclear effects considered here and often neglected in the Monte-Carlo 
CIhI simulations, the absolute events rates are subject to important modifications. 
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1 Introduction 

Neutrino physics is among the hottest topics of particle physics with the recent indications 
in favor of neutrino oscillations. After the solar neutrino deficit, the apparent anomaly in the 
ratio of muon to electron atmospheric neutrinos R^/e = {Ny^ + Np^)/{Ny^ + NpJ observed 
by (Super-)Kamiokande [^ |[, 1MB [Q], Soudan-2 0] and the asymmetry in the zenithal 
distributions of the /i — type events in Super-Kamiokande have given a strong support to 
the oscillation hypothesis: z/^ — > v^ where v^ = Vt (^-e. active-active transition) or z/^ = Vg 
[i.e. active-sterile transition). The solution of the atmospheric neutrinos anomaly in terms 
of z/^ — >■ Ve oscillations has been excluded by the Chooz collaboration . 

A number of atmospheric neutrinos experiments use large underground water Cerenkov 
detectors. In these experiments only "one Cerenkov ring" (1 C.R.) events are retained for 
the analysis. These events are usually assumed to be produced by quasi-elastic charged 
current interactions in which a charged lepton is emitted above Cerenkov threshold and 
leads to one Cerenkov ring. The nucleon which is ejected from the nucleus is in general 
below threshold and therefore does not produce another ring. The region of energy transfer 
in processes involving atmospheric neutrinos of ~ 1 GeV extends from the quasi-elastic peak 
to the Delta resonance region. The evaluation of the nuclear responses in the latter region 
usually relies on the assumption that the Delta decays into a pion and a nucleon (this is 
the case for example in ref. where the authors use a relativistic model a la Walecka to 
compute the nuclear response functions). The pion leading to an additional Cerenkov ring. 



this charged current event belongs to the two Cerenkov rings (2 C.R.) class and is rejected 
by the experimental cuts. Thus theoretical calculations are often limited to the quasi-elastic 
peak which is treated in Fermi gas models or with more elaborate treatments taking into 
account the shell structure of the oxygen nucleus and RPA type correlations . 

However the nuclear dynamics is far more complex than this simple picture. Indeed the 
pion in the nucleus is a quasi-particle with a broad width and can decay for instance into 
a particle-hole excitation. Therefore the decay of a Delta in the nuclear medium can lead 
to a nucleon and a particle-hole state. In such a process, two nucleons are ejected from the 
nucleus, none of them producing a Cerenkov ring, and the event belongs to the 1 C.R. class. 
Furthermore {2p-2h) states may also be directly excited in the nucleus without excitation of 
the Delta resonance. This process also results in the emission of two nucleons and the event 
belongs to the 1 C.R. class. Following these arguments, we perform a full calculation of the 
neutrino-oxygen cross sections beyond the quasi-elastic assumption, with the identification 
of the possible final states. This procedure leads to a complete evaluation of the 1 C.R. 
events yields in the atmospheric neutrinos experiments and its impact in the description of 
the retained neutrino events in the detectors has to be investigated. 

The starting point of this calculation is the inclusive charged current cross section for 
the reaction z/; {vi) +^^ O — > l' (/+) + X, 
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where Gp is the weak coupling constant, 6c the Cabbibo angle, k and k' the initial and final 
lepton momenta, q^ = k^ — k' = {uj, q) the four momentum transferred to the nucleus, 9 the 
scattering angle. Ma (M) the Delta (nucleon) mass. The plus (minus) sign in eq. (|ID stands 
for the neutrino (antineutrino) case. In a provisional approximation, to be lifted after, we 
have neglected in eq. (P the lepton masses and we have kept the leading terms in the devel- 
opment of the hadronic current in p/M, where p denotes the initial nucleon momentum. The 
electric, magnetic and axial form factors are taken in the standard dipole parameterization 
with the following normalizations: Ge{^) = 1-0, Gm{^) = 4.71 and G'a(O) = 1.25. We have 
introduced the inclusive isospin (Rr), spin-isospin longitudinal [R„^(^i^) and spin-isospin 
transverse {RariT)) nuclear responses functions (the longitudinal and transverse character 
of these last two responses refers to the direction of the spin operator with respect to the 
direction of the transferred momentum) : 

^r = E (^1 E Or(j) e^^-^^ |0)(n| E O^J{k) e^9-^M0)* 6{u -E^ + E^) (2) 

n i=l fc=l 

where the operators have the following forms: 

OaU) = T-f. {(^j-Q) rf, {{(Tj xq)xq) rf , 



for a = r, aT{L), ar{T), and 

0^{j) = {S,.q)T^, {{S,xq)xq)T^, 

for a = aT{L), aT{T). In the above expressions, the superscript P {P = N or A) denotes the 
type of the Particle-hole excitations {Nucleon-hole or Delta-hole) induced by the operator 
O^. The operators S and T are the usual 1/2 to 3/2 transition operators in the spin and 
isospin space (for instance see |^). In this work we neglect the small quadrupole transition 



connecting the nucleon to the Delta through the pure isospin operator, therefore the isospin 
response just involves nucleon-hole excitations. Note that we have assumed the existence of 
a scaling law between the nucleon and Delta magnetic and axial form factors [p!0[] : 

Ml^M — (JAI^A — J I h 

where /* (/) is the it N A (tt N N) coupling constant. For a matter of convenience, we have 
incorporated the scaling factor /*// = 2.2 into the responses. 

2 Formalism 

The evaluation of the nuclear responses is performed within the model developed by Delorme 
and Guichon for the interpretation of the ('^He,t) charge exchange experiments [|1T], |I2|. In 
this model the polarization propagators Il^{uj,q,q') without nuclear correlations are evalu- 
ated in a semi-classical approximation to properly take into account the finite size effects. 
This implies the use of a local Fermi momentum kp{r) which is calculated by the means of 
an experimental nuclear density: A;^(r) = (3/2 7r^p(r))^/^. Note that this procedure differs 
a little from the pure semi-classical one but it has been found to give better results for the 
71 — nucleus reactions. The "bare" polarization propagators 11° (in the following "bare" will 
mean that the nuclear correlations are switched off) are then used as an input to exactly 
solve the RPA equations in the ring approximation, as we will develop in the following. In 



ref. [0, |I^ the authors gave satisfactory fits to the set of the experimental data. This model 
was also confronted to the pion-nuclei experimental results [jl3| and the agreement obtained 
for the total and elastic cross sections was remarkably good. 

As mentioned above, the first step of the calculation is the evaluation of the bare polariza- 
tion propagators. A crucial ingredient of the model is the Delta resonance width modified by 



the nuclear effects. We adopt the parameterization of ref. |]T3 where the Delta width is split 



into the contributions of different decay channels: the "quasi-elastic" channel, A — > n N, 
modified by the Pauli blocking of the nucleon and the distortion of the pion, the two-body 
{2p-2h) and three-body {3p-3h) absorption channels. This parameterization leads to a good 
description of pion-nuclear reactions. At resonance we find a Delta width around 130 MeV, 
a value rather close to the free case. This value reflects the importance of the two- and 
three-body absorption channels which are large enough to counteract the effect of the Pauli 
blocking and lead to this overall enhancement of the Delta width. Note furthermore that 
at resonance the "quasi-elastic" channel modified by the medium effects, is almost equal 
to the free "quasi-elastic" one. The model of Delorme and Guichon also accounts for the 
{2p-2h) excitations which are not reducible to a modified Delta width. The evaluation of 

^For a pure quantum approach in the low energy part of the nuclear response, applied in the context of 
terrestrial neutrinos experiments, see ref. [|l Bl 



such processes is performed by extrapolating the calculations of two-body pion absorption 
at threshold given in ref . [^ . We have limited ourselves to the imaginary part of these two- 
body polarization propagators, the comparison with experimental data such as pion-nucleus 
scattering or (e, e') scattering giving satisfactory results to that order of approximation. By 
construction, the bare polarization propagator Il^{uj, q, q') is the sum of the following partial 
components: 

1. A^A^ quasi-elastic (the standard Lindhard function), 

2. A^A^ {2p-2h), 

3. A^A and 3'. AA^ {2p-2h), 

4. AA (ttN), 

5. AA {2p-2h), 

6. AA i3p-3h), 

where the notation A^ (A) stands for Nucleon-hole (Delta-hole) states as previously. The 
Feynman graphs corresponding to this partial polarization propagators are displayed on 
fig. ([1|) with the following conventions: the wiggled lines represent the external probe, the 
full lines correspond to the propagation of a nucleon (or a hole), the double lines to the 
propagation of a Delta, the dashed lines to an effective interaction between nucleons and/or 
Deltas. Finally the dotted lines indicate which intermediate state has to be placed on-shell 
to obtain the desired partial nuclear response. Note that in the case of [np-nh) polarization 
propagators the number of graphs is large and we just give one example in the figure. 




Figure 1: Feynman graphs of the partial polarization propagators: NN quasi-elastic (1), 
NN (2p-2h) (2), NA (2p-2h) (3), AN (2p-2h) (3'), A A (nN) (4), A A (2p-2h) (5), A A 
(3p-3h) (6). The conventions for the various lines drawings are given in the text. 



The bare responses are then given by the standard relations: 



R\k){uj,q) = —\m{Yi\k){uj,q,q)), 

n 



(3) 



with the obvious sum rule: 

R\u,q) = —lm{U\u;,q,q)) = J2 R\k){u;,q), (4) 

^ fc=i 

where n^ denotes the number of partial reaction channels ( n^ = 7 in our model). 

Following the method detailed in ref. [|ri|, |r^ we include the effects of nuclear correlations 
by exactly solving the RPA equations in the ring approximation. For instance the inclusive 
RPA polarization propagators Il{uj, q, q') are solution of the generic equation: 

n = n° + n° 'i/n (5) 

where V denotes the effective interaction between particle-hole excitations and n°((x', q, q') 
the inclusive bare polarization propagator calculated previously and used here as an input. 
In the spin-isospin channel the RPA equations couple the L, T and the iV, A components 
of the polarization propagators. For the effective interaction relevant in the isospin and 
spin-isospin channels, we use the standard vr + p + 5 — function parameterization: 

VnN = (/' + V^ + Vp + Vg.) T1.T2 

Vna = (K + V, + Vg,) T,.Tl 
Van = (K + V, + Vg,) T,.T2 
Vaa = (K + V, + Vg,) T,.Tl (6) 

where in the A^A^ case, for example (the A^A, AA^ and A A cases are obtained with the 
appropriate replacements a — > S): 



K = F^ (^ ^ t) (^1-^ CTa-g 

„2 

Vp = Fl (— -) <Ti X g 0-2 X g 

V,, = Fl g' CT,.CT, (7) 

In the preceding equations, -F^(g) and Fp{q) are the standard pion-nucleon and rho-nucleon 
form factors. The values we adopt for the relevant parameters can be found in ref. [jT^. In 
particular we take the "common" Landau- Migdal parameters: 

/' = 0.6, g'j^^ = 0.7, g'^^ = g'^^ = 0.5, g'^^ = 0.5 

The inclusive RPA responses functions are deduced from the corresponding inclusive RPA 
polarization propagators by the usual relation: 

R{LU,q) = --lm{U{u;,q,q)) (8) 

TV 

We perform the calculation of the partial RPA responses as follows. Starting from the RPA 
equation (j^), we write the imaginary part of Il{uj, q, q') in the following form: 

im(n) = |i + nv"|2im(n°) + ini^imV" (9) 



with Im(n°) = J2kLi^^^^{k))- This sum rule gives the different contributions to the 
inclusive RPA response functions. The first terms in eq. (^ are reminiscent of the bare case. 
Indeed we recognize the bare partial response functions (apart from the trivial — tt factor) 
corrected by a factor involving the inclusive RPA polarization propagator and the effective 
interaction. The partial RPA response functions, defined by: 
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(10) 



are represented by the graphs (a) to (d) on fig. (|^), where the hatched rings correspond to 
the inclusive polarization propagator solution of eq. (^), the non hatched rings to the bare 
partial polarization propagators (the dotted line means that we take the imaginary part of 
these propagators) and the dashed lines to the effective interaction. It is easy to recover on 
these graphs the different terms of the development of eq. (p^Of). Note that in the RPA case, 
a PP' reaction channel (P, P' = N, A) gets contributions from every QQ' configurations 
(A^A^, A^A, AA^, AA). The last term in eq. (|^) corresponds to the "coherent" response 
function: 

R,,h{uJ,q) = --\n\HmV (11) 



n 



It is absent of the response spectrum when the effective interaction is switched off. In the 
domain of energy considered here the sole contribution to this channel comes from the pion 
exchange. This process corresponds to the emission of a pion on its mass-shell, the nucleus 
remaining in its ground state. It is represented by the graph (e) on fig. (^ where the dashed 
line stands for the exchange of a pion. The implications of these partial reaction channels 
will be discussed in the following sections. 




Figure 2: Graphic representation of the partial RPA response functions, (a)-(d): incoher- 
ent partial response functions, (e): coherent partial response function. The hatched rings 
correspond to the inclusive RPA polarization propagator, the free rings to the partial hare po- 
larization propagators, the dashed line to the effective interaction and the dotted line indicates 
which intermediate state is placed on-shell. 



3 Cross sections 

The next step is the calculation of the neutrino-oxygen cross section. The doubly differential 
cross section is given in a first approximation by eq. (|l|). Our final calculation relies on 
a more complete expression, which we will briefly describe in the following, but the main 
features remain unchanged. First it is essential to note that the neutrino-nucleus reaction 
is strongly dominated by the transverse spin-isospin channel. This is clear from eq. (|l]). 
Indeed the terms multiplying the transverse responses, depending on the axial and magnetic 
form factors, have a much larger magnitude than for the longitudinal case. Furthermore 
the A^A^ quasi-elastic spin-isospin longitudinal response is totally suppressed in the cross 
section. This suppression arises from an exact cancellation at the top of the quasi-elastic 
peak between the various terms entering the contraction of the leptonic and the hadronic 
tensors. This suppression is only partial in the Delta resonance or in the "dip" region (the 
region intermediate between the quasi-elastic and the Delta peaks). This result is in con- 



tradiction with the study of ref. [16| where the relative weights of the transverse to the 
longitudinal responses were assumed to be 2:1. Note however that the suppression of the 
NN longitudinal response is no more exact when one considers the complete expression of 
the doubly differential cross section, which includes the terms involving the charged lepton 
mass (in fact we consider only the muon mass) and the terms up to order (p/M)^ in the 
reduction of the hadronic current. The contributions of the NN longitudinal response are 
then of order [mi/My, where rui denotes the mass of the charged lepton, and of order {p/M). 
These corrections are rather weak. As another source of corrections we have also considered 
the renormalization of the axial charge by the mesonic exchange currents, because the sup- 
pression of the A^A^ longitudinal response involves the time component of the axial current. 



Following the parameterization of ref. |17| we make the replacement gA — ^ 5'a(1 + ^) in the 
time component of the axial current. The contribution of the A^A^ longitudinal response is 
then of order P. Even with the relatively high value S ~ 0.5 suggested in ref. [|18|, the 



contribution of the A^A^ longitudinal response remains weak. The same conclusion holds for 
the A^A and AA longitudinal responses which are widely dominated by the corresponding 
transverse ones. Note that in the antineutrino-nucleus reactions the weight of the transverse 
channel is somewhat reduced because of the change in sign in the interference term. However 
even in this case, the transverse responses correspond to 75 % of the total, the remaining 
arising essentially from the A^A^ pure isospin response. 

We will now investigate the implications of these global features on the simply differen- 
tial cross section da/dk', which is obtained from the doubly differential cross section by a 
numerical integration over the solid angle. The great interest of our method is the separation 
of the inclusive cross section on partial contributions. This separation is simply achieved by 
the replacement, in the expression of the cross section, of the inclusive response functions 
with the "exclusive" ones, calculated in the previous section. The results are shown on fig. 
(j^) which displays the simply differential cross section versus the energy transfer, fixing for 
the sake of illustration a neutrino energy of 700 MeV. 

The inclusive cross section is given by the thick curve. It gets its main contribution from 
the A^A^ quasi-elastic channel (thin full line) which peaks at relatively low energy transfer. 
For the sake of comparison we have shown the contribution of the A^A^ quasi-elastic channel 
without RPA (thin long dashed hne). We observe that the cross section is reduced and 
hardened in the RPA case. This result is in full agreement with that of ref. |]^. The 
shift in strength reflects the dominance of the transverse response. Indeed the Landau- 
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Figure 3: Differential charged current z/g — ^^ O interactions cross-section versus the energy 
transfer. The thick curve represents the inclusive cross-section. The following exclusive con- 
tributions to the inclusive cross-section are displayed: NN (quasi-elastic) (full thin curve), 
NN (2p-2h) (short dotted curve), NA + AA^ (2p-2h) (short dot-dashed curve), AA (tcN) 
(short dashed curve), AA (2p-2h) (long dot-dashed curve) and AA (3p-3h) (long dotted 
curve). Also shown is the "hare" NN quasi-elastic cross-section (long dashed curve). 



Migdal interaction is repulsive for all values of the transfer and the p — exchange piece is not 
attractive enough in the domain of energy considered here to counteract this feature. This 
repulsive effective interaction hardens and reduces the transverse response functions. For 
instance this conclusion no longer holds in the longitudinal channel where the vr — exchange 
is attractive enough to overcome the Landau-Migdal interaction and to create a collective 
mode (the so-called pionic branch [5,6]). But, as we mentioned previously, the suppression 
of the longitudinal channel makes the neutrino a poor probe of this pionic mode. 

The effect of the RPA correlations are less strong in the others channels and are not shown 
here. The AA (vrA^) channel (short dashed curve on the figure) arises at high energy transfer 
{uj ~ 450 MeV). This is in good agreement with the result of the relativistic calculations 
given in ref. [^] where the Delta was taken into account as a free resonance and where the 
RPA correlations did not include Delta-hole configurations. The agreement between the two 
calculations is not surprising. Indeed we use relativistic kinematics for the evaluation of the 
polarization propagators as in ref. [5,6] and we include the terms up to second order in the 
(p/M) reduction of the hadronic current as mentioned above. Furthermore the RPA effects 
on the transverse response in the Delta region, that is at high transfer, are somewhat reduced 
in the cross section by the form factors and the differences between RPA and bare transverse 
response functions are very weak at these values of the transfer. Thus our calculations show 



that a free Delta resonance gives a rather good approximation of the A A (ttN) (or "quasi- 
elastic") channel. This result corroborates the fact that the "quasi-elastic" Delta width in 
the nuclear medium is close to the its free value, the Pauli blocking being cancelled by the 
other mechanisms taken into account. 

The most interesting feature of the cross section is the importance of the {np-nh) channels. 
The kinematics of the neutrino- nucleus reaction tends to favor the NN {2p-2h) channel 
(short dotted curve) which peaks at low energy transfer. However the A^A -|- AA^ {2p-2h) 
channel (short dot-dashed curve) gives a rather large contribution to the inclusive cross 
section and has an extended spectrum in the " dip" region. Its importance has been pointed 
out in the (e, e') scattering where it is necessary to reproduce the experimental data in the 
"dip" region (for example see ref. |19[)- Finally note that the AA {2p-2h) (long dot-dashed 



curve) and {3p-3h) (long dotted curve) spectra extend over a wide range of energy transfer, 
while the AA (ttN) channel is concentrated in the so-called Delta peak. They give a little 
contribution to the inclusive cross section (in particular the {3p-3h) channel is rather weak) 
but the extension of their spectra will have important consequences in the specific events 
yields. 

The results obtained in this section show that the inclusive neutrino-oxygen cross section 
is strongly modified with respect to the free NN quasi-elastic case, which is quite often the 
sole channel entering into the calculations. In particular we have seen the occurrence of 
large contributions from the two- and three-body channels. The main effect of the nuclear 
correlations is the hardening of the A^A^ quasi-elastic channel. They have rather low impact 
on the others reaction channels and therefore could be legitimately neglected. 

4 Events yields 

In this section we compute numerically the neutrino-oxygen events yields for a fixed charged 
lepton momentum: 



Y{u^ + u^){k') = J2 dE {^.AE)^,{E,k')^,SE)^,{E,k') 



(12) 



where E is the neutrino energy, (^y ($p) the incoming neutrinos (antineutrinos) flux and 
da/dk' [da/dk') the neutrino-oxygen (antineutrino-oxygen) cross section computed in the 
previous section. We use the fluxes of Bartol pO| in our calculations for the sake of compar- 
ison with the results of ref. |^ . The main feature of these fluxes is their sharp decrease with 
increasing neutrino energy. Note that several theoretical attempts have been undertaken 
to compute these atmospheric neutrinos fluxes. The sources of possible differences between 
three models have been analyzed in ref. [^. The predictions of these models on the flavor 



ratio agree at a ~ 5 % degree of accuracy. Anyway the divergences in the absolute fluxes 
remain rather large (~ 20 %). Furthermore new measurements on the primary cosmic rays 
fluxes could lead to some modifications with respect to the present situation. The cumu- 
lated uncertainties on the neutrino fiuxes and on the neutrino-oxygen cross sections could 
lead to modifications in the experimental analysis, even if they remain unlikely to explain 
the atmospheric neutrinos anomaly. To perform an analysis of the events yields, we need 
to classify the partial reaction channels according to the number of Cerenkov ring(s) they 
produce. This classification has been elaborated within a few rough assumptions. First we 
consider that every nucleon ejected from the nucleus remains under Cerenkov threshold. In 



'O charged current 1 C-ring interaction rates 



water, the threshold kinetic energy for a particle of mass m is E r^ 0.5 m. For a nucleon, 
produced through Delta decay or {2p-2h) mechanisms, the assumption is fairly good. On 
the opposite, we assume that every pion which escapes the nuclear medium produces a 
Cerenkov ring. The threshold energy for a pion being ~ 70 MeV, this assumption is believed 
to be reliable. Then the partial reaction channels leading to one Cerenkov ring, in charged 
current interactions, are the A^A^ quasi-elastic one, which is usually taken into account in 
the simulations, and the {np-nh, n>l) type channels (both NN, NA, AN and AA). The 
remaining reaction channels, AA (vrA^) and coherent pion production, are supposed to lead 
to at least two Cerenkov rings. The results for the 1 C.R. events yields, which are relevant 
in the atmospheric neutrino experiments, are shown on fig. (^ for incident z/^ and z/^, the 
full curves corresponding to the total 1 C.R. events yields and the dashed curves to the sole 
A^A^ quasi-elastic 1 C.R. events yields. We give the results of the calculations without (thin 
curves) and with (thick curves) RPA. 

First we observe that the RPA tends 
to reduce the events yields. This is 
not hard to understand. Indeed the 
RPA tends to harden the cross sec- 
tions, i.e. to push the strength towards 
higher energies. But the fluxes decrease 
with increasing energies and therefore 
the higher energies are disfavored. This 
reduction affects mostly the A^A^ quasi- 
elastic channel in accordance with the 
result obtained for reaction cross sec- 
tions. The maximal reduction factor is 
of the order of 10 %. A more interest- 
ing feature is the strong enhancement 
of the absolute events yields imphed by 
the (np-nh) channels. At the maximum 
value of the yields, the enhancement of 
the total yield with respect to the A^A^ 
quasi-elastic one is around 30 %. This 
result reflects the main features of the 
cross sections. Furthermore one must be 
aware that this result is a lower limit of 
the true enhancement. Indeed we know 
that pions can be re-absorbed in the nu- 
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Figure 4: 

yields versus the muon momentum, 
curves correspond to the total 1 C.R. events yields 
with (thick curve) and without (thin curve) RPA. 
The dashed curves correspond to the NN quasi- 
elastic 1 C.R. events yields with (thick curve) and 
without (thin curve) RPA. 

cleus. Therefore the events produced in the AA (vrA^) channel can also lead to one Cerenkov 
ring if the pion does not escape from the nucleus. Thus we can conclude that the RPA 1 C.R. 
events yields induced by charged current interactions is globally enhanced with respect to the 
A^A^ quasi-elastic 1 C.R. events yields without RPA. The difference between the two calcula- 
tions could be responsible for the small discrepancy between the experimental and simulated 
events distributions in Super-Kamiokande 0. But we need complementary informations to 
ensure this conclusion. Indeed it is hard to establish the enhancement factor firmly. Our 
present analysis leads to a maximum enhancement factor of the order of ~ 20 %. We have 



^We thank Y. Declais for attracting our attention on this point. 



already mentioned that the problem of pion absorption, which is not yet considered in our 
calculations, could still enhance this factor. We should also be aware that there exists some 
misidentification problems which could have a more or less large effect. One of the misiden- 
tification source, pointed out by the authors of ref. 0, is the "coherent" pion production. 
Their analysis is based on the assumption that the forward peaked angular distribution of 
the coherent pions entails the coherent pions to be emitted with a small angle with respect 
to the charged lepton direction. This could mimic a "shower" which could be interpreted as 
an e — type event, whatever may be the flavor of the incoming neutrino. Our calculations 
show that this coherent channel brings a tiny contribution (less than 2 % of the total z/g + Ve 
events yield) which makes it irrelevant in the atmospheric neutrino anomaly. The suppres- 
sion of this channel is understandable. Part of it is due to the nuclear form factors effects 
as discussed for example in [^. In addition the coherent response manifests itself mainly in 
the longitudinal spin-isospin channel and we have seen that this channel is suppressed with 
respect to the transverse one in the neutrino-nucleus reactions. The coherent pions should 
not be a problem. 

The case of the neutral currents is less clear. In charged current interactions, pions 
("coherent" or not) lead to, at least, 2 C.R. events and are excluded from the analysis. But 
in neutral currents interactions they lead to 1 C.R. events, because the scattered neutrino 
does not produce any ring. We have computed the neutral current events yields in each 
reaction channel, the few differences with respect to the charged current case being easily 
included in the formalism. The problem arising then is the classification of these vr — like 
1 C.R. events. Indeed in absence of indication on the experimental vr/lepton discrimination 
efficiency in the water Cerenkov experiments, it is not possible to draw firm conclusions on 
the role played by the neutral currents. This problem has to be investigated further (for 
more details on the results of the calculations, see [p^). 

Finally we study the evolution of the flavor ratio with the charged lepton momentum: 

R^,e{k') = Y{v^ + P^)/F(z/, + z^e), (13) 

where Y denotes the events yields defined by eq. (PTSj). We compare the total 1 C.R. events 
yields fiavor ratio with the A^A^ quasi-elastic 1 C.R. events yields fiavor ratio. The result 
is shown in table (|l]) where the ratio of ratios has been calculated for four relevant lepton 
momentum. There is almost no modification between the two situations. This conclusion 
strengthens the usual assumption that uncertainties due to nuclear effects cancel when one 
considers ratios of events rates. Here the maximum effect on the flavor ratio is less than 
10%. 



fc' (MeV/c) 


i?^/e(iViVg.e.)/i?^/e(Tota/) 


100 


1.060 


150 


1.040 


250 


0.999 


400 


1.001 



Table 1: Comparison of the total and NN quasi-elastic 1 C.R. events yields ratios for four 
lepton momentum. 



We conclude this work by mentioning the problem of pion emission in neutrino-oxygen 
interactions. On one side we have shown that the cross sections of the AA {2p-2h) and {3p- 



3h) partial channels, which do not not lead to pion emission {non pionic channels), extend 
over a broad region in transfer energy, while the pionic channel AA (ttN) is peaked at high 
transfer energy (see fig. @). On the other side the neutrino flux lowers the weight of the 
high energies and favors the low energy components of the spectrum. Then the pionic AA 
channel will be more suppressed by the incident neutrino flux than the non pionic one. This 
result is shown on fig. (^ where the total AA events yield (full thick curve) is split into its 



three contributions: (vrA^) (full thin curve) 
the case of /i — type events. 

The main result is that the fraction 
of the non pionic channels over the pi- 
onic one is around 50 %. This result 
remains valid for every values of the 
lepton momentum. Thus the (np-nh) 
excitations play an important role in 
the events yields although their reaction 
cross sections are relatively low. Finally 
we would like to point out that some 
pion production mechanisms, which do 
not reduce to a simple response func- 
tion, are still absent of our formalism. 
For example we have omitted in the pro- 
duction through the vector current, the 
Kroll-Ruderman and the pion-in-fiight 
terms which play an important role. Im- 
provements to our present calculations 
are in progress. Nevertheless these lim- 
itations of our present calculations do 
not alter the need of including the ef- 
fects of the partial AA (np-nh) reaction 
channels to avoid an overestimation of 
the number of pions effectively produced 
in the neutrino-oxygen interactions. 



-2h (dashed curve), 3p-3h (dotted curve) in 
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Figure 5: Contributions to the total AA events 
yield (full thick curve) of the partial channels: 
(ttN ) (full thin curve), (2p-2h) (dashed line), (3p- 
3h) (dotted line). 



5 Conclusion 



In this work we have studied the effects of nuclear correlations on the charged current 
neutrino-oxygen cross sections and events yields in specific exclusive reaction channels. We 
have shown that besides the quasi-elastic channel the {np-nh, n=2,3) excitations also lead to 
one Cerenkov ring events, which are retained for the analysis of the experiments using large 
underground water Cerenkov detectors. The enhancement in the one Cerenkov ring events 
yields is large and could still be increased when some processes, such as pion absorption in 
nuclei or neutral currents events, are taken into account. It is therefore important to take 
these nuclear effects into account to perform a calculation of absolute events yields. We have 
also shown that the flavor ratio R^/e is not significantly altered. We have applied our model 
to others neutrino-nucleus reactions. In particular we have studied the case of iron which 
is the target-nucleus in the neutrino experiments using calorimeters. The conclusion on the 



cross sections are the same than the one presented here in the case of oxygen. However 
such experiments measure more detailed observables than the water Cerenkov detectors, hke 
the energy and momentum spectra of the particles in the final state. The description of 
these experiments requires the extension of our model. The present work already shows the 
necessity of taking into account nuclear correlations involving multi-nucleon excitations. 
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